Raman scattering from a superconductivity-induced bound state in MgB 2 
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It is shown that the sharp peak in the E2 g Raman spectrum of superconducting MgBi is due to 
a bound state caused by the electron-phonon coupling. Our theory explains why this peak appears 
only in the spectra with i?2 9 symmetry and only in the a but not it bands. The properties of the 
bound state and the Raman spectrum are investigated, also in the presence of impurity scattering. 
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Electronic Raman scattering in superconductors 
probes the density and, to some extent, the momen- 
tum dependence of quasi-particle excitations across the 
gap. As a result the spectra are usually dominated by 
broad pair breaking peaks reflecting the magnitude and 
anisotropy of the gap as well as scattering processes of the 
quasi-particles. Experimental spectra in the two-gap su- 
perconductor MgBi show the typical features of a dirty 
s-wave superconductor [1,2]. A new and surprising fea- 
ture in the Ei g spectrum is a sharp peak near the larger of 
the two gaps which in our opinion cannot be explained 
as a pair breaking peak because of the large impurity 
scattering rate of quasi-particles in the present MgB 2 
samples. Below we will argue that this peak is the ana- 
logue of the magnetic resonance peak observed in high-T c 
cuprates [3] where the electron-phonon interaction takes 
over the role of the Heisenberg interaction and strongly 
scatters the excitations across the gap. 

Neglecting the momenta of the incident and scattered 
photons the differential Raman cross section can be writ- 
ten as 



= -iL(i + „(„))*» 



in 



e is the charge of an electron, n(tu) the Bose distribu- 
tion function, u> the difference in frequency between in- 
cident and scattered light, and x" the imaginary part of 
a retarded Green's function x a t zero momentum. The 
corresponding Matsubara function is 
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where (3 is the inverse temperature 1/T, uj n — 2nnT 
are bosonic Matusbara frequencies, T T the time ordering 
operator, and p an effective density operator with mo- 
mentum zero. In the following we will be interested in 
spectra related to the E-ig symmetry of D§hi the point 
group of MgB%. p has then the form 
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cj. n (k) and Co- n (k) are creation and annihilation opera- 
tors for electrons with momentum k, spin direction a, 
band label n, and energy e n (k). MgB^ has only one 
Raman-active q = phonon, and this phonon has Ei g 
symmetry. Its frequency is denoted by f2, its creation 
and annihilation operators by V s and b, respectively, and 
the corresponding element of the phononic Raman ten- 
sor by R. Without loss of generality we assume that 
this phonon transforms as the first basis vector of the 
two-dimensional E 2g representation in accordance with 
Eq.(4). 

For the evaluation of x we use the Hamiltonian H = 
H + H' with 

H Q = e„(k)ct n (k)c CTn (k) + ]T ^(q)(&](q)&,(q) + \) 
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A„ is the gap parameter for s-wave superconductivity in 
the band n, £lj (q) the frequency of a general phonon with 
momentum q and branch label j, g the coupling constant 
for intraband electron-phonon scattering and V a ran- 
dom potential for intraband impurity scattering. Inter- 
band phonon scattering can be neglected in H' because 
only zero momentum transfers occur in the approxima- 
tion used below. 

In a first step we perform an infinite summation over 
bubble diagrams by introducing the irreducible Green's 
function x- X contains all diagrams to x which cannot be 
decomposed into two parts by cutting one phonon line. 
The average over impurities yields impuritiy lines with 
similar properties as phonon lines. However, bubble di- 
agrams connected by impurity lines are not possible in 
this case so that the above definition of irreducibility is 
appropriate. Analytically, one obtains 
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X = Xii + (X12 + R)D(R + X21), 
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The omitted frequency and momentum arguments in the 
Green's functions in Eqs.(7) and (8) are iuj n and 0, re- 
spectively, xn denotes the irreducible Green's function 
associated with the two vertices 7„(k). Similarly, the ver- 
tices in the functions X12 and %22 are 7n(k) and <?„(k0), 
and two times g n (k0), respectively. The free phonon 
progator is given by — 2Q/(uj 2 + Q 2 ). 

A sensible approximation for the evaluation of \ is the 
ladder approximation plus the corresponding self-energy 
corrections where the interaction lines are due to phonons 
or impurities. Only that part of the interaction can con- 
tribute in the ladder diagrams which transforms in the 
same way as the vertices 7 and g. This means in our case 
that only the Ei g component of the phonon-mediated in- 
teraction, which usually is considered to be neglegible in 
a s-wave superconductor, would enter. We therefore will 
evaluate \ orL ly m the presence of random impurities us- 
ing the Born approximation and the dirty limit for each 
band. Assuming that the vertices and the interaction 
can be evaluated right on the Fermi surface, expanding 
7 and g in terms of Fermi surface harmonics [4] of E 2g 
symmetry, <]>^(k), normalized as in Eq.(l) of Rcf. [4], 
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and assuming that the interaction is diagonal in L we 
obtain 

XiA^n) =J2a^a^N^ x (iu; n ,A n ,r- 1 ). (11) 

Ln 

(n) 

Np is the density of states at the Fermi surface for 
one spin direction due to the band n and 1/t„ an ef- 
fective scattering rate. After an analytic continuation 
iu> n — > u> + ir\ the imaginary part of the Green's function 
X, which is independent of L, becomes [5] 
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is the theta function and F, E, and II are complete 
elliptical integrals of the first, second, and third kinds, re- 
spectively. The existence of only two gaps in the experi- 
mental spectra [1,2,6] as well as theoretical arguments [7] 
suggest that the dirty limit applies even within the two- 
band complex of a and tt bands (denoted by the index p 
in the following) and that the interband impurity scat- 
tering between a and tt bands is neglegible. As a result 
A„ and r" 1 can be considered to be the same within the 
manifold of a or tt bands. Introducing then the effective 
couplings 
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we obtain 
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Choosing for the first function $^™' > (k) the properly nor- 
malized function ~ 7n(k), the summations over L col- 
lapse to one term L = 1 if either i or j is equal to 1. In 
the clean limit 1/t — > Eq.(12) reduces to the analytical 
formula Eq. (16) of Ref.( [8]). 

Using the tight-binding fit to the band structure of 
Ref. [9] , eV as energy and the lattice constant a as length 



units, we find Njp 



N 



(4) 



0.049, 



0.109, N ( p ] = 0.101, 



0.104, where n = 1, 2 denote the light and heavy 



a bands and n = 3,4 the tt bands. In the case of tt 
bands we have scaled the published tight-binding pa- 
rameters slightly in order to reproduce the correct densi- 



ties. Furthermore, we obtain = 0.101, ol( } = 0.261 



„( 2 ) 
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2.531, a{ i] = 2.376, and thus X^' = 0.008 



and A^ = 1.238. The rather small values for a\ 1! , 
and A^ reflect the fact that the a bands are 
rather isotropic in the ab-plane near their narrow cylin- 
drical Fermi surfaces and thus cannot contribute much 
in the E 2g channel. As a consequence we may safely put 
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X^i = 0. The a contribution to x arises 
then solely from the coupling of the phonon to the light 
and its spectral features are determined by the dimen- 
sionless electron-phonon coupling constant A = ^A^ 
and the scattering rate I/tv. 

We will first consider the a contribution to x and drop 
the index a everywhere in order to simplify the notation. 
Fig.l shows the real (dashed lines) and imaginary (solid 
lines) part of % for the gap 2A = 110cm -1 and two scat- 
tering rates 1/r = and 1/t = 200cm~ 1 . In the clean 



2 



case the imaginary part x" exhibits a square root sin- 
gularity if to approaches 2A from above. Impurity scat- 
tering transforms this singularity into a step at 2A with 
height 27tAt, and produces for sufficiently strong scat- 
tering a very broad minimum near 1/r. The real part 
X' shows in the clean case a square root singularity be- 
low 2 A and finite and positive values above 2 A. Taking 
impurity scattering into account x' becomes much more 
symmetric around 2A compared to the clean case, 



lu£ = n 2 + {i 2 \x'(u b ,A,i/ T ) 
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FIG. 1. Real (dashed lines) and imaginary (solid lines) part 
of the susceptibility \ f° r two impurity scattering rates r -1 
and parameters appropriate for MgB^- 

exhibiting there a logarithmic singularity and negative 
values over a large region above 2 A. 
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FIG. 2. Energy oj(,/2A (solid curves) and spectral weight 
Zf, (dashed curves) of the bound state for two impurity scat- 
tering rates r -1 as a function of the electron-phonon coupling 
constant A. 



The curves in Fig. 1 suggest that the phonon Green's 
function D will develop a bound state inside the gap [10], 
and that this will occur for all values for A and 1 jr. The 
frequency of the bound state, u)b, is determined by the 
vanishing of the denominator of D, i.e., by the equation, 



Expanding the denominator of D around iOb one finds for 
the spectral weight Z b 



2Q 



(19) 



Fig. 2 shows u>b and Zb as a function of A for the two 
scattering rates 1/r = and 1/r = 200cm -1 . u>b ap- 
proaches zero at Xcdw = 0.50 and 0.47, respectively. 
This means that for A > Xcdw the superconducting 
state is unstable against the formation of a charge den- 
sity wave with E2 g symmetry. For A < Xcdw ^>b and 
Zb increase and decrease rapidly with decreasing A ap- 
proaching their limiting values 1 and in an exponential 
(for 1/r ^0) or powerlawlike (for 1/r = 0) manner. 

One peculiar feature of D is that any background imag- 
inary part in its self-energy, for instance due to ir elec- 
trons, will be diminuished by the factor 1/Zb near 2 A. 
This means that the bound state will not be broadened 
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FIG. 3. Susceptibility x" °f G electrons at T = for the 
scattering rate r" 1 = 200cm -1 . 

and disappear in x" f° r small A's, even in the presence 
of strong impurity scattering, but instead will sharpen 
and loose weight when ujb approaches 1. To make the 
situation more realistic one should allow for broaden- 
ing effects due to inhomogenieties etc. which may be 
taken into account by folding x' with a Gaussian with a 
width S [2,5]. Fig. 3 shows x" f° r three different A and 
a width of 5 = 10cm -1 . The employed large scattering 
rate r -1 = 200cm -1 has wiped out completely the usual 
pair breaking peak due to the square root divergence of 
x" at 2 A. On the other hand, the electron-phonon cou- 
pling accumulates spectral weight near 2A and produces 
a pronounced bound state inside the gap at larger cou- 
plings. 

Figs. 2 and 3 allow to estimate a realistic value for A. 
Identifying LUb with the observed sharp peak in the Ei g 
spectrum we have uib = 104 ± 1 cm -1 whereas a re- 
cent tunneling experiment [6] gave for the gap 2A = 
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114 ± 6 cm . This suggest that lui, is different from 
2 A, i.e., the sharp peak should not be identified with the 
gap, and the ratio is 0.91 ± 0.06. Fig. 2 indicates 

then that A must be smaller than 0.3. Comparing the in- 
tensities of the bound state and the phonon line in Fig. 3 
with the experimental curve [1] one finds that A must be 
near the interval between 0.15 and 0.20. These values are 
substantially smaller than the value 0.38 obtained from 
band structure calculations [11]. 

Additional evidence for the importance of the electron- 
phonon coupling in the E 2g spectrum comes from the ex- 
perimental result that the two A ig and the E lg spectra 
do not show any pronounced peak near 2A. We explain 
this by the fact that, according to group theory, MgB 2 
has noq = phonons with such symmetries so that no 
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FIG. 4. Experimental (solid line) and theoretical (dashed 
line) E2 g Raman cross section. 

bound states can be formed. Further support for the 
present theory comes from the absence of a peak near 
the gap of about 47 cm -1 in the it bands. LDA cal- 
culations show that the coupling of n electrons to the 
E2 g phonon is by a factor 3 or 4 smaller than for a 
electrons [11]. We find that the bound state structure 
near the 7r gap is invisible for such a coupling which ex- 
plains its absence in the experimental spectra. Using 
X[l> = A&> = 0, AM/4? = 1/3, X^/R 2 = V400, 
A = 0.16, £1 = 620cm 1 , S = 5cm 1 , a 7r gap of 43cm 1 
and a somewhat reduced a gap of 96.5 cm -1 , account- 
ing for the finite temperature of the experimental data, 
the resulting Raman cross section is shown as the dashed 
line in Fig. 4, together with the experimental E 2g spec- 
trum [1]. The dashed line reproduces well the main fea- 
tures of the experimental curve, especially at small fre- 
quencies. The quantitative discrepancy in the phonon 



region suggests that only part of the phonon broaden- 
ing is caused by the electron-phonon coupling and that 
anharmonicity and deviations from the assumed con- 
stant density of states may play a role. The theoretical 
superconductivity-induced hardening of the phonon fre- 
quency in Fig. 4 is 5 cm -1 and thus somewhat smaller 
than the experimental value [1] of 7 — 10 cm -1 . 

In conclusion, we have shown that the E 2g spectrum in 
superconducting MgB 2 can be understood as a superpo- 
sition of a phonon line coupled strongly to a electrons cre- 
ating hereby a bound state in the gap, and a background 
due to rather uncorrelated tt electrons. This means that 
MgB 2 is to our knowledge the first s-wave superconduc- 
tor where a bound state in the superconducting gap due 
to residual interactions has been observed and identified. 
The obtained electron-phonon coupling constant A <~ 0.2 
for a electrons is only half of the band structure value, a 
discrepancy, which presently is not well understood. 

The author thanks O. Dolgov, J. Kortus and I. Mazin 
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